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We develop a numerical algorithm to calculate the degrees of irreducible multiparty correlations 
for an arbitrary multiparty quantum state, which is efficient for any quantum state of up to five 
qubits. We demonstrate the power of the algorithm by the explicit calculations of the degrees of 
irreducible multiparty correlations in the 4-qubit GHZ state, the Smolin state, and the 5-qubit W 
state. This development takes a crucial step towards practical applications of irreducible multiparty 
correlations in real quantum many-body physics. 
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Introduction. — In modern physics, as is well known, 
mean field theory is not sufficient to describe the physics 
in a strongly correlated many-body system [H, which im- 
plies that there exist rich correlation structures in its 
quantum state. Therefore, how to characterize multi- 
party correlations in a multipartite quantum state be- 
comes a fundamental problem in many-body physics. 
Traditional physical method is to introduce the corre- 
lation functions to describe the correlations in a many- 
body system. 

The extensive researches on characterizing entangle- 
ment in quantum information science [2] shed different 
light on the problem. On one hand, correlation functions 
are not invariant under local unitary transformations, 
which implies that they can be regarded only as correla- 
tion witnesses but not as legitimate correlation measures 
0, 0|- On the other hand, the information based view- 
point can be instructive in characterizing correlations in 
a multipartite quantum state. 

In the information viewpoint, the degree of the to- 
tal correlation in a multipartite quantum system is 
equal to the difference between the sum of von Neumann 
entropies of all the subsystems and the von Neumann 
entropy of the whole system. There are two different 
schemes to classify the total correlation: one is to distin- 
guish the total correlation into quantum correlation and 
classical correlation 0, 0] , the other is to divide the total 
correlation into pairwise correlation, triplewise correla- 
tion, etc.. 

The concept of irreducible n-party correlation in an 
n-partite quantum state was first proposed in Ref. Q. 
We generalized it to irreducible m-party (2 < m < N) 
correlation in an n-partite state, and proposed that all 
the irreducible m-party correlations construct a classifi- 
cation of the total correlation @, 0|. It is worthy to note 
that, in classical information community, the irreducible 
m-party correlations in a joint probability distribution 
of n classical random variables were investigated in Ref. 

13, EH- 

The degree of irreducible multiparty correlations in a 
multipartite quantum state, like many important quan- 
tities in quantum information science, such as the mea- 



sure of entanglement [12 | and the capacity of a quantum 
channel ljj 14 1, are defined as an optimization problem, 
which makes their calculations become extremely diffi- 
cult. These computational difficulties almost prevent any 
practical application of these measures in a real physical 
problem. Therefore it is of great significance to develop 
an efficient algorithm to calculate them for a general mul- 
tiparty quantum state. 

In Ref. Q, we proposed a continuity approach that 
reduces the calculations of irreducible multiparty corre- 
lations in a multiparty quantum state without maximal 
rank to the calculations of irreducible multiparty correla- 
tions in a series of multiparty quantum states with max- 
imal rank. Although theorem 1 in Ref. Q tells us the 
form of a maximal rank state without higher order ir- 
reducible multiparty correlations, however, this theorem 
does not solve the problem on the calculations of irre- 
ducible multiparty correlations for a general multiparty 
quantum state with maximal rank. This is why we solve 
the calculations only for some specific classes of states 
in Ref. ||. In other words, we have not a systematic 
method to calculate the irreducible multiparty correla- 
tions for a state with maximal rank. 

In this Letter, we develop an efficient systematic nu- 
merical algorithm on the calculations of the degrees of 
irreducible multiparty correlations for a general multi- 
partite quantum state. The advantages of our algorithm 
is that it is independent of initial values of variables, and 
we find it is efficient for an arbitrary quantum state of up 
to five qubits. To the best of our knowledge, it is for the 
first time that we have the capacity to deal with the de- 
tailed analysis of the correlations in a general multiparty 
state of up to five qubits. 

Notations and Definitions. — The Hilbert space of 
an n-partite quantum system is denoted by = 
n™ = i®W (l) , where [n] is the set {1,2, •••,«}, and 
is the Hilbert space of party i whose dimension is 
di. The inner product of two operators and B^> 
in the Hilbert space is defined as (A®\B®) = 

^-Tt(A^B^) [13]. The prefactor £ is introduced 
to satisfy the normalization condition = 1, 
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where jW is the identity operator in the Hilbert space 
Hi. Thus we can introduce an orthonormal Hermitian 
operator basis {Oi,,a, € {0, 1, — 1}}. In par- 

ticular, we take Oq to be the identity operator /W. 
Any operator can be expanded in this basis as 

Ail) = J2a z ° ( ^(Oa}\A^). Further more, the opera- 
tor A^l in the n-party Hilbert space can be expanded 
as AM = £ a(n) O^O^AM), where a(n) is the set 

{ai, aa, ■ • • , a n }, and O^"^ is an abbreviated notation for 

n™=i ^o?- If the operator AM is the Hamiltonian of an 
n-party system, the terms to describe m-party interac- 
tions (1 < m < n) satisfy the condition iVo(a(n)) = n — m 
with JV (a(n)) = <W 

Without loss of generality, we consider an n-party 
quantum state pM with maximal rank, which can be ex- 
panded as 

p iBi =E°!?n,(oa,iP iB] >. (i) 

a(n) 

Because the state pM is positive definite, we can define 
InpM uniquely as a Hermitian operator. Then we can 
apply the above expansion to In pM to obtain 

^^E ^ ^^}- (2) 
a(n) 

The condition TrpM = 1 implies that the coefficient 
^o(ri) I m ) can ^ e determined by the other coefficients 
(Oi" ] n) | InpW). Here 0(n) is the set a(n) with a, = for 
if [n], and a(n) is the same as a(n) except 0(n). Com- 
pared with the expansion JT]), the obvious advantage of 
the expansion ([2]) is that it ensures the positivity of pM 
automatically. Further more, a one-to-one map between 
the state pM with maximal rank and the set of real co- 
efficients { (0~™'„) | In pM } } can be built. The existence 
of such a one-to-one map is an essential element in our 
numerical algorithm. 

To make use of the expansion |(2]) , we adopt the equiv- 
alent definitions of the degrees of irreducible multiparty 
correlations in a multiparty quantum state given in Ref. 
[l| but not the original definitions proposed in Ref. 
If we adopt the original definition, then the optimization 
is made under the expansion {1} , which makes the opti- 
mization almost impossible because of the constraint of 
semipositivity of a density matrix. In Ref. 0, we give 
the definitions on the degrees of irreducible multiparty 
correlations for a three-qubit system. Now the defini- 
tions are generalized for a general multipartite quantum 
state with a finite dimensional Hilbert space as follows. 

We first define the set of the n-party states without 
more-than-m-party irreducible correlations as 

B m = {<r [n] 1(0^1 InaW) = 0,VJV o ( a (n)) < n - m}. 

(3) 



Next we find the state in the set B m that is least distin- 
guishable with the state pM 

pN^arg min S(p^\\a^), (4) 

where the quantum relative entropy S(p\\p') = 

Tr(p(lnp — In//)) for two quantum states p and p' in 
the same Hilbert space. Then the degree of irreducible 
m-party correlation is defined as 

^(pHj^^HpNj. (5) 
In addition, the degree of the total correlation is defined 

by 

CT{pW)=S{pW\\f}p). (6) 

Using the same arguments given in Ref. [jj, we can show 
that C T (pM) = E: =2 C m {pW) = £?=i S(pU)-S(pW) 
with the von Neaumann entropy S(p) = — Tr(plnp) for 
a quantum state p. 

Numerical Algorithm. — In the above optimization 

In] 

problem, it is an essential task to find out the state p m for 
a given state pM, It is possible to directly solve Eq. (f4]) 

In] 

to obtain the state p m . However, it is doutful whether 
the solution we find is a local minimum or a global min- 
imum. Fortunately, the optimization problem |(4]) can be 
transformed into the following system of nonlinear equa- 
tions: 

<Oi"L)I lll < ™> =0 ' V7V (a(n))<n-m, (7) 
(O l $ n) \p% ] ) = (0 [ $ n) \p [n] ), VN ( a (n))>m. (8) 

In Ref. we proved that there exists a unique real 
solution of {(OK L | hi Pm )} satisfying the above system 
of equations for a three-qubit system. This result is also 
valid for a general multipartite quantum state with a fi- 
nite dimensional Hilbert space. Here we neglect the proof 
because it is a simple generalization for the three-qubit 
case. Thus we have two different ways to use the system 
of equations (0 HJ. On one hand, we can use them to 
verify whether the solution of the optimization problem 
(0]) is correct. On the other hand, we can directly use the 
optimization method to solve them to obtain the states 
Pm- In our present numerical algorithm, we adopt the 
latter method in application of the system of equations 
(HI [3]). We want to emphasize that Eq. is used to 
represent a multipartite quantum state in our algorithm. 

For an optimization problem, one of the key skills is to 
choose a proper initial value. Here we adopt a continuity 
approach to choose a proper initial value for any n-partite 
quantum state pM . We consider a seriers of states 

Tin] 

P W (Po)=Po^ R + (l-Po)p [ " 1 . (9) 
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We take p = 1 - ^ with k <E {0, • • ■ , N}, where N is 
a large positive integer. Obviously, (fc — 0) = , 
pM(fe = N) = pM and (0'™ ] Jin p^ffc = 0)) = 0. We 
take the values of { (Og/L | hi p™ (A;)) } as the initial values 
of {(0™ n) | ln^ 1 (k + 1))} for k = 0, 1, • • • , N - 1. 

The basic idea under the above approach is based 
on the continuity principle, more precisely, the state 
+ 1) is very similar to the state p^(fc), so the 
values of {(OW J lupin (k + 1))} is also near the values 

{(Og(L| lnpm (fc))}. The practice of our computations 
shows that our selection of initial values makes the algo- 
rithm become efficient. The cost of the algorithm is that 
we calculate the degrees of irreducible multiparty corre- 
lations for a series of states p'"'(po) instead of a single 
state p [n] . 

An obvious advantage is that the choose of initial val- 
ues in our algorithm is independent of the state p^ . In 
other words, our algorithm makes the computation of a 
general multiparty state become efficient. In my personal 
computer, it is efficient for any state up to five qubits. To 
the best of my knowledge, it is the best results on mul- 
tiparty correlations in a multiparty state we obtained so 
far. 

Numerical results. — We will demonstrate the power of 
our numerical algorithm by explicitly giving the results 
on irreducible multiparty correlations for some typical 



multiparty states: the 4-qubit GHZ state [17(1, * ne f° ur - 
qubit Smolin state [la ], and the 5-qubit W state 

The first state we consider is the n-qubit GHZ state 
|GHZ n ) = ^ (nr=i ®|0}* + Uti Th e degrees of 

irreducible multiparty correlations on the 4-qubit GHZ 
state are given in Figure 1. The total correlation in the 
state is 4 bits, and it is classified into 3 bits of irreducible 
two-qubit correlation and 1 bit of irreducible four-qubit 
correlation. These results are the same as those given in 
Ref. @L and they are consistent with the conclusion in 
Ref. (20, ED]. 

The second state we consider is the 4-qubit Smolin 
state, whose density matrix is simply given by p^mo = 

Te( ili] +Ui=i +n<Li We find that there exists 

2 bits of correlation in the state, and they are irreducible 
4-qubit correlations, which is shown in Figure 2. From 
the density matrix of the Smolin state, we know that it 
is also a generalized stabilizer state defined in Ref. 0. 
In this sense, the numerical results also verify the results 
in Ref. [|. 

The third state we consider is the 5-qubit W state 
| W5) = Ya=i 11^=1 Our numerical results 

show that only irreducible two-qubit correlations exist in 
the W state, which support the conclusion in Ref. ([22]). 

In the range of our numerical results, we find that 
the degree of the total correlation Ct is a non-increasing 
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Figure 1: The degrees of irreducible multiparty 
correlations for the 4-qubit GHZ state. 




Figure 2: The degrees of irreducible multiparty 
correlations for the four-qubit Smolin state 
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Figure 3: The degrees of irredcible multiparty 
correlations of the 5-qubit W state. 
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function of po, however, the degree of irreducible m-party 
correlation can increase with increasing po (see, for exam- 
ple, Figure 3). Actually we can prove that Ct{p^ (po)) 
is a non-increasing function of po for any n-party state 
as follows. We can imagine that every subsystems of 
the n-partite quantum system pass through a depolarized 
channel [23| , then the quantum state evolves accord- 
ing to Eq. lO in the direction of increasing po. In the 
process, only local operations act on the state p™, and 
the degree of total correlation does not increase under lo- 
cal operations, therefore Ct(p^(po)) is a non-increasing 
function of po for any n-party state p^ . 

In addition, the fact C m (m = 3,4,5) is not a non- 
increasing for a 5-qubit W state gives another example 
to support one of the main results in Ref. [1]: local 
operations can transform lower order correlations into 
higher order correlations. 

Discussions and summary. — The calculations of the 
degrees of irreducible multiparty correlations for an arbi- 
trary multiparty quantum state are challenging because 
they are defined as the constraint optimization prob- 
lems over all the multiparty quantum states in the whole 
Hilbert space. In this Letter, we develop an efficient nu- 
merical method to calculate the degrees of irreducible 
multiparty correlations for any multipartite quantum 
state, which is based on the following two key elements. 

One key element in our algorithm is that we adopt 
the expansion of a multipartite state in the exponential 
form ([2]). First, it ensures the positivity of the state 
automatically. Second, although the independent vari- 
bles {(OW, | lnpl™])} can take the limit to infinity, the 

state is always well defined because of the constraint 
Trpl n l = 1. In this sense, the state without maximal 
rank are naturally contained in this expansion if the co- 
efficients can limit to infinity. This makes our algorithm 
effective for arbitrary multiparty states. 

The other key element is related to the selection of the 
initial values of variables, more precisely, the formula iJH). 
It makes our algorithm independent on the initial values 
of variables, and greatly enhances the efficiency of our 
algorithm. 

In summary, we present an efficient numerical algo- 
rithm on the calculations of the degrees of irreducible 
multiparty correlations in a multiparty quantum state. 
Our algorithm is valid for arbitrary quantum states up 
to five qubits in my personal computer, and it is a univer- 
sal algorithm whose efficiency does not depend strongly 
on the multipartite quantum state. We demonstrate the 
power of our algorithm by explicitly giving the results for 
the 4-qubit GHZ state, the Smolin state, and the 5-qubit 
W state, which are consistent with previous results. We 



hope that our development of this algorithm will provide 
a powerful tool to analyze the correlation distributions 
in a multipartite quantum state, and thus takes a cru- 
cial step towards the practical applications of irreducible 
multiparty correlations in real quantum many-body sys- 
tems. 

The author thank Dr. S. Yang and Prof. CP. Sun 
for stimulating discussions. This work is supported by 
NSF of China under Grant No. 10775176, and NKBRSF 
of China under Grants No. 2006CB921206 and No. 
2006AA06Z104. 



[1] F. Alet, A.M. Walczak, and M.P.A. Fisher, Physica A 

369, 122 (2006). 
[2] R. Horodecki, P. Horodecki, M. Horodecki, and K. 

Horodecki, Rev. Mod. Phys. 81, 865 (2009). 
[3] L. Henderson and V. Vedral, J. Phys. A 34, 6899 (2001). 
[4] D.L. Zhou, B. Zeng, Z. Xu, and L. You, Phys. Rev. A 

74, 052110 (2006). 
[5] S. Watanabe, IBM Joural of Research and Development 

4, 66 (1960). 

[6] B. Groisman, S. Popescu, and A. Winter, Phys. Rev. A 

72, 032317 (2005). 
[7] N. Linden, S. Popescu, and W.K. Wootters, Phys. Rev. 

Lett. 89, 207901 (2002). 
[8] D.L. Zhou, Phys. Rev. Lett. 101, 180505 (2008). 
[9] D.L. Zhou, Phys. Rev. A 80, 022113 (2009). 
[10] S. Amari, IEEE Trans. Inf. Theory 47, 1701 (2001). 
[11] E. Schneidman, S. Still, M.J. Berry II, and W. Bialek, 

Phys. Rev. Lett. 91, 238701 (2003). 
[12] C.H. Bennett, D.P. DiVincenzo, J. A. Smolin, and W.K. 

Wootters, Phys. Rev. A 54, 3824 (1996). 
[13] A.S. Holevo, IEEE Trans. Inf. Theory 44, 269 (1998). 
[14] B. Schumacher and M.D. Westmoreland, Phys. Rev. A 
56, 131 (1997). 

[15] R.A. Horn and C.R. Johnson, Matrix analysis (Cam- 
bridge University Press, 1985). 

[16] V. Vedral, Rev. Mod. Phys. 74, 197 (2002). 

[17] D.M. Greenberger, M.A. Home, and A. Zeilinger, Bell's 
theorem, Quantum Theory, and Conceptions of the Uni- 
verse (Kluwer Academics, Dordrecht, 1989). 

[18] J. A. Smolin, Phys. Rev. A 63, 032306 (2001). 

[19] W. Diir, G. Vidal, and J.I. Cirac, Phys. Rev. A 62, 
062314 (2000). 

[20] S.N. Walck and D.W. Lyons, Phys. Rev. Lett. 100, 
050501 (2008). 

[21] N. Linden and W.K. Wootters, Phys. Rev. Lett. 89, 
277906 (2002). 

[22] P. Parashar and S. Rana, Phys. Rev. A 80, 012319 
(2009). 

[23] M.A. Nielsen and I.L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
New York, 2000). 



